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ABSTRACT 

The nearby low-luminosity active galactic nucleus (LLAGN) NGC 4258 has a weak radio continuum 
component at the galactic center. We investigate its radio spectral properties on the basis of our 
new observations using the Nobeyama Millimeter Array at 100 GHz and archival data from the 
Very Large Array (VLA) at 1.7-43 GHz and the James Clerk Maxwell telescope at 347 GHz. The 
NGC 4258 nuclear component exhibits (1) an intra- month variable and complicated spectral feature 
at 5-22 GHz and (2) a slightly inverted spectrum at 5-100 GHz (a ~ 0.3; F v oc v a ) in time- 
averaged flux densities, which are also apparent in the closest LLAGN M81. These similarities between 
NGC 4258 and M81 in radio spectral natures in addition to previously known core shift in their AU- 
scale jet structures produce evidence that the same mechanism drives their nuclei. We interpret the 
observed spectral property as the superposition of emission spectra originating at different locations 
with frequency-dependent opacity along the nuclear jet. Quantitative differences between NGC 4258 
and M81 in terms of jet/counter jet ratio, radio loudness, and degree of core shift can be consistently 
understood by fairly relativistic speeds (F > 3) of jet and their quite different inclinations. The picture 
established from the two closest LLAGNs is useful for understanding the physical origin of unresolved 
and flat /inverted spectrum radio cores that are prevalently found in LLAGNs, including Sgr A*, with 
starved supermassive black holes in the present-day universe. 

Subject headings: galaxies: active — galaxies: individual(NGC 4258) — galaxies:jets — galaxies: 
Seyfert — radio continuum: galaxies — Submillimeter: galaxies 



1. INTRODUCTION 

Most of the supermassive black hole activity in the 
present-day universe is primarily in the lower end of 
the luminosity function of active galactic nuclei (AGNs) 
in the form of low-luminosity AGNs (LLAGNs; for re- 
views, see Ho et al. 1997; Ho 2008) with bolometric 
luminosities Lboi ^ 10 42 erg s _1 down to nearly qui- 
escent systems such as Sgr A* at our Galactic cen- 
ter. Most LLAGNs are extremely sub-Eddington sys- 
tems (Lt>oi/-i'Edd < 10~ 3 , where £ecM is the Eddington 
luminosity). The low luminosities and spectral energy 
distributions (SEDs) of LLAGNs (Ho 1999; Eracleous 
et al. 2010; Younes et al. 2012) have been modeled by 
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the combination of (i) an inner advection-dominated ac- 
cretion flow (ADAF; Narayan & Yi 1994), which is a 
radiatively inefficient, optically thin, and geometrically 
thick accretion flow at low accretion rates; (ii) an outer 
truncated disk (Quataert et al. 1999) of a standard opti- 
cally thick, geometrically thin accretion disk (Shakura & 
Sunyaev 1973); and (hi) a jet (Falcke & Biermann 1995; 
Yuan et al. 2005), as proposed by several studies (e.g., 
Nemmen et al. 2006, 2011; Yu et al. 2011, and references 
therein) . 

The ADAF model predicts an SED including a "sub- 
millimeter bump" that results from the synchrotron 
emission of thermal hot (~ 10 9 -10 10 K) electrons. A 
highly inverted spectrum (a = 0.4; Mahadevan 1997, in 
S F a where a is the spectral index and S„ is the flux 
density at frequency v) is expected in lower-frequency ra- 
dio bands. However, observed spectra in LLAGNs tend 
to be flat or slightly inverted (a ~ —0.2 to +0.2) in the 
central components (Nagar et al. 2001; Doi et al. 2005b, 
2011) as well as the spectrum of Sgr A* at < 10 GHz (Fal- 
cke et al. 1998; An et al. 2005). Furthermore, only the 
thermal process in the ADAF model cannot provide ob- 
served radio luminosities adequately in many cases (Fab- 
biano et al. 2003; Doi et al. 2005a; Wu & Cao 2005; Wu 
et al. 2007). Thus, these are the reason why the jet 
component is needed in SED modeling for LLAGNs. It 
is also argued that the entire SED of LLAGNs is prac- 
tically jet-dominated by synchrotron emission and up- 
scattered inverse Compton component rather than by the 
accretion flow or disk (in particular at very low luminosi- 
ties Lx/^Edd % 10~ 6 ; Yuan & Cui 2005), according to 
studies for individual sources (NGC 4258; Yuan et al. 
2002a, M81; Markoff et al. 2008) and many samples in 
the statistical sense (Merloni et al. 2003; Falcke et al. 
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2004; Kording et al. 2006; Yuan et al. 2009; de Gasperin 
et al. 2011; Plotkin et al. 2012). Thus, the jet component 
is presumably essential in the energetics at the low state 
of AGNs. 

Very-long-baseline interferometry (VLBI) images at 
milli-arcsecond (mas) resolutions have sometimes re- 
vealed elongated radio structures similar to sub-parsec 
scale jets in LLAGNs (Falcke et al. 2000; Nagar et al. 
2005). With the exception of a handful of well-known 
nearby radio galaxies and Seyfert galaxies with jets ex- 
tending to kpc scales (Wrobel 1991; Ho & Ulvestad 2001), 
most of the power is concentrated in an unresolved core 
of < 10 3 Rs-10 4 Rs (Rs is the Schwarzschild radius) with 
high brightness temperatures of < 10 6 -10 11 K (Ulves- 
tad & Ho 2001; Anderson et al. 2004; Filho et al. 2004; 
Krips et al. 2007). Radio flux variability observed in sev- 
eral LLAGNs in total flux on typical timescales of a few 
days also indicates < lOOORs in size (Anderson & Ulves- 
tad 2005). Thus, the radio emitting origin has been still 
unknown for most LLAGNs from an observational point 
of view. One of the rare cases in which a nuclear radio 
component has been spatially resolved is M81 (Bieten- 
holz et al. 1996), which is the nearest (3.6 Mpc; Freed- 
man et al. 2001) LLAGN in the VLBI-observed sample 
of Nagar et al. (2005). M81 contains a type 1.5 Seyfert 
nucleus (Ho et al. 1997) with a 2-10 keV X-ray lumi- 
nosity of <~ 1.5-4 x 10 40 erg s _1 (Ishisaki et al. 1996; 
Markoff et al. 2008) and a black hole with a mass of 
7 x 10 7 M Q (Devereux et al. 2003). The radio nucleus 
consists of a bright core plus weak one-sided jet-like elon- 
gation with a scale of <~ 1 mas (3 x 10 3 Rs) at 1-43 GHz in 
VLBI images (Bietenholz et al. 2004; Ros & Perez- Torres 
2012). Observed radio properties are characterized by as 
follows. (1) Core shift: the position of radio brightness 
peak appears shifted depending on observing frequency 
at 1.7 8.4 GHz; the line of sight toward the putative 
black hole is opaque at these radio frequencies (Bieten- 
holz et al. 2004; Marti- Vidal et al. 2011). Therefore, al- 
most all radio fluxes should be attributed to AU-scale jets 
rather than an accretion flow. (2) Flux variability: M81 
is the most well-studied LLAGN in total flux as well; 
a bright radio nucleus (~ 70-400 mJy) exhibits rapid 
and large amplitude intraday fluctuation (Ho et al. 1999; 
Sakamoto et al. 2001) and on the timescale of several 
weeks and years (Bietenholz et al. 2000; Markoff et al. 
2008; Marti- Vidal et al. 2011). (3) Inverted spectrum: 
a centimeter-to-submillimeter (cm-to-submm) spectrum 
is slightly inverted on time average (a ~ 0.3; Reuter & 
Lesch 1996; Doi et al. 2011) with a possible turnover fre- 
quency at 100-230 GHz (Schodel et al. 2007). (4) Spec- 
tral variability: a complexly variable spectral profile was 
observed at meter to submm wavelengths (235 MHz- 
345 GHz) in quasi-simultaneous observations over six 
months (Markoff et al. 2008). Thus, the radio properties 
(2)-(4) at arcsec resolutions are ascribable to nuclear jets 
whose structure is spatially resolved as (1). For this rea- 
son, M81 may be a Rosetta Stone for understanding the 
origin of unresolved radio cores in other LLAGNs and 
Sgr A*. 

NGC 4258 contains a type 1.9 Seyfert nucleus (Ho 
et al. 1997) with a 2-10 keV X-ray luminosity of (4.2- 
17.4) xlO 40 erg s _1 (Makishima et al. 1994; Fruscione 
et al. 2005) and a black hole with a mass of 3.9 x 10 7 M 
(Herrnstein et al. 1999). NGC 4258 is the second near- 



est (7.2 Mpc; Herrnstein et al. 1999) LLAGN in the 
VLBI-observed sample of Nagar et al. (2005). Thereby, 
NGC 4258 in addition to M81 constitute rare cases 
in which a nuclear radio structure is spatially resolved 
(Miyoshi et al. 1995). VLBI imaging at - 1.5 GHz 
showed a two-sided nuclear jet extending in the north- 
south direction (Cecil et al. 2000) that appears phys- 
ically connected with the kpc-scale jet-like radio struc- 
tures ("anomalous arms"; e.g., van der Kruit et al. 1972). 
The nuclear radio structure has been thoroughly inves- 
tigated at 22 GHz through VLBI observations of water 
masers (Herrnstein et al. 2005; Argon et al. 2007). A 
two-sided nuclear jet with a scale of <~ 1 mas (9 x 10 3 Rs) 
upwells from the dynamical center of a putative disk of 
Keplerian rotating water masers in a disk with a nearly 
edge-on orientation. VLBI images provided a 3cr upper 
limit of 220 [iJy on 22 GHz at the dynamical center; most 
of the nuclear radio continuum flux comes from the jet, 
which is located significantly off the dynamical center 
(Herrnstein et al. 1998). The jet exhibits a flux variation 
of <~ 100% on the timescales of a few weeks or longer 
(Herrnstein et al. 1997, see also Appendix A). Thus, 
the observed properties in terms of innermost jet struc- 
ture and flux variability of the nuclear component are 
very similar to those of M81 ((1) and (2)). However, the 
radio spectral properties are known in much less detail 
(Appendix A): (3) a slightly inverted spectrum? (4) a 
variable spectral profile? To establish the archetype of 
LLAGNs by using M81 and NGC 4258, the radio spectral 
properties of NGC 4258 nucleus should be investigated 
and compared with those of M81 at cm-to-submm wave- 
lengths. 

In this paper, we present the radio continuum spec- 
trum and its variability for the NGC 4258 nucleus in 
the cm-to-submm bands at arcsecond resolutions. This 
study also includes imaging of extended components, 
such as extended synchrotron jets and dust emission in a 
host galaxy, for the purpose of estimation of contamina- 
tion. The observations and data reductions are described 
in Section 2. The results are presented in Section 3. The 
physical origins of observed radio properties and compar- 
isons with M81 are discussed in Section 4. Finally, we 
summarize in Section 5. For the distance of 7.2 Mpc to 
NGC 4258, 1" corresponds to 35 pc. 

2. OBSERVATIONS AND DATA REDUCTIONS 

For the purpose of revealing radio spectral variability 
and its average nature for this variable radio nucleus in 
NGC 4258, we obtained multi-epoch data. We carried 
out new observations of 15 epochs using the Nobeyama 
Millimeter Array (NMA) at the Nobeyama Radio Ob- 
servatory (NRO) at - 100 GHz (Section 2.1). We re- 
trieved the 1.7-43 GHz archival data of the Very Large 
Array (VLA) at the National Radio Astronomy Obser- 
vatory (NRAO) for 13 epochs in total (Section 2.2). The 
list of these observations is presented in Table 1. We 
reduced the 347 GHz archival data of the James Clerk 
Maxwell telescope (JCMT) on Mauna Kea (Section 2.3), 
which was, as a result, devoted to estimation of dust 
contamination in photometry toward the active nucleus 
because the nuclear region was dominated by dust emis- 
sion at this frequency (Section 3.1). 

2.1. NMA Observations 
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TABLE 1 

Data List of VLA and NMA Observations 



Date 


Project 


Array 


V 






P.A. 


a 




code 


config. 


(GHz) 


(") 


(") 


O 


(mjy beam -1 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


07 Doc 1996 


AH594A 


VLA-A 


22.460 


0.10 


0.09 


34 


0.12 








14.940 


0.15 


0.14 


-8 


0.10 








8.460 


0.29 


0.23 


-12 


0.04 








4.860 


0.47 


0.43 


-18 


0.06 


22 Dec 1996 


AH594B 


VLA-A 


22.460 


0.11 


0.09 


28 


0.11 








14.940 


0.15 


0.13 


-7 


0.11 








8.460 


0.28 


0.23 


-14 


0.05 








4.860 


0.47 


0.43 


-21 


0.06 


29 Doc 1996 


AH594C 


VLA-A 


22.460 


0.10 


0.09 


22 


0.13 








14.940 


0.15 


0.14 


-2 


0.11 








8.460 


0.28 


0.24 


-11 


0.06 








4.860 


0.46 


0.44 


-10 


0.07 


31 Dec 1996 


AH594D 


VLA-A 


22.460 


0.10 


0.09 


29 


0.17 








14.940 


0.15 


0.14 


-6 


0.15 








8.460 


0.24 


0.20 


-15 


0.05 








4.860 


0.47 


0.43 


-19 


0.06 


04 Jan 1997 


BR0043 


VLA-A 


15.365 


0.16 


0.14 


-87 


0.11 








8.415 


0.70 


0.27 


-66 


0.07 








4.985 


0.92 


0.44 


-69 


0.07 


05 Jan 1997 


AH594E 


VLA-A 


22.460 


0.10 


0.09 


28 


0.18 








14.940 


0.16 


0.14 


-6 


0.15 








8.460 


0.29 


0.23 


-13 


0.05 








4.860 


0.47 


0.43 


-20 


0.07 


07 Jan 1997 


BG0062 


VLA-A 


1.665 


1.23 


1.41 


83 


0.04 








22.460 


0.17 


0.11 


-69 


0.23 


06 Mar 1997 


BM0056 a 


VLA-B 


22.265 


0.35 


0.31 


-66 


0.10 








8.460 


1.88 


0.83 


75 


0.06 








4.860 


3.54 


1.43 


69 


0.12 


19 Mar 1997 


AH594F 


VLA-B 


22.460 


0.43 


0.32 


82 


0.20 








8.460 


1.42 


0.90 


85 


0.12 








4.860 


2.24 


1.47 


-83 


0.14 


28 Feb 1998 


AG0527 


VLA-A 


43.340 


0.05 


0.04 


-13 


0.41 








22.460 


0.11 


0.09 


86 


0.33 


05 Sep 1998 


BM0112 a 


VLA-B 


22.265 


0.40 


0.33 


87 


0.16 








14.940 


1.97 


0.45 


56 


0.16 








8.460 


2.41 


0.81 


62 


0.07 








4.860 


3.74 


1.41 


64 


0.10 


21 Dec 2000 


AN0097 


VLA-A 


43.340 


0.09 


0.04 


76 


0.51 








8.460 


0.42 


0.24 


79 


0.11 


21 Dec 2001 




NMA-D 


100.777 


7.5 


6.2 


-72 


0.89 


21 Mar 2002 




NMA-C 


100.777 


8.8 


4.7 


-36 


3.00 


22 Dec 2003 




NMA-C 


100.777 


4.3 


3.5 


-2 


0.94 


14 Jan 2004 




NMA-AB 


100.777 








1.96 


31 Mar 2004 




NMA-D 


100.777 


6.6 


5.7 


-76 


0.93 


01 Apr 2004 




NMA-D 


100.777 








2.76 


05 Apr 2005 




NMA-C 


100.777 


3.8 


3.0 


-33 


0.84 


08 Apr 2005 




NMA-C 


100.777 


4.4 


3.6 


-30 


1.49 


13 May 2005 




NMA-D 


95.729 


7.5 


6.2 


-88 


0.87 


14 May 2005 




NMA-D 


95.729 


9.5 


7.9 


-70 


1.93 


15 May 2005 




NMA-D 


95.729 


9.9 


6.8 


76 


0.94 


25 Mar 2006 




NMA-D 


100.777 








3.12 


29 Mar 2006 




NMA-D 


100.777 








3.12 


31 Mar 2006 




NMA-D 


100.777 








2.37 


01 Apr 2006 




NMA-D 


100.777 








2.86 



Note. — Column 1: observation date; Column 2: project code; Column 3: telescope and array configuration; Column 4: observing frequency; 
Columns 5—7: FWHMs in major and minor axes for synthesized beam and position angle of the major axis; Column 8: image rms noise. 

a VLA participated in VLBI observation as a phased-up element, 
and the VLBI results have been published (Argon et al. 2007). 



We observed the nuclear regions at <~ 100 GHz 
(A3 mm) by using the NMA. Array configurations were 
AB, C, and D with maximum baseline lengths of 351, 
163, and 82 m, respectively. The visibility data were ob- 
tained with double-sided receiving systems. We use only 
the upper sideband at a center frequency of 100.777 GHz 



with a bandwidth of 1 GHz, which was practically line 
free; the emission lines such as HCN(J =1-0) and 
HCO + (J =1- 0) were in the lower sideband at a cen- 
ter frequency of 88.777 GHz. For the data acquired 
in 2005 May, we use two sidebands centered at 89.729 
and 101.729 GHz, both of which were practically line 
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TABLE 2 

Flux Densities of Radio Continuum Emission at the Nucleus of NGC 4258 



Date Si. 7 S5 58.4 S15 S22 S43 S100 S347 q 5 

(mjy) (mjy) (mjy) (mjy) (mjy) (mjy) (mjy) (mjy) 

W (2) (3) (4) (5) (6) (7) (8) (9) (10) 

07 Dec 1996 1.8 ± 0.1 1.9 ± 0.1 2.2 ± 0.2 2.4 ± 0.3 0.21 ± 0.04 

22 Dec 1996 2.0 ± 0.1 2.2 ± 0.1 2.2 ± 0.2 2.2 ± 0.2 0.06 ± 0.06 

29 Dec 1996 1.9 ± 0.1 2.3 ± 0.1 2.7 ± 0.2 2.5 ± 0.3 0.21 ± 0.07 

31 Dec 1996 1.9 ± 0.1 2.3 ± 0.1 2.7 ± 0.3 3.1 ± 0.4 0.32 ± 0.01 

04 Jan 1997 1.8 ± 0.2 2.0 ± 0.1 3.5 ± 0.2 0.64 ± 0.23 

05 Jan 1997 2.3 ± 0.1 2.3 ± 0.1 2.6 ± 0.3 3.6 ± 0.4 0.19 ± 0.12 

07 Jan 1997 2.0 ± 0.1 b 3.7 ± 0.4 

06 Mar 1997 3.1 ± 0.3 2.7 ± 0.2 5.4 ± 0.3 0.62 ± 0.19 
19 Mar 1997 2.7 ± 0.3 2.8 ± 0.3 3.2 ± 0.4 0.10 ± 0.01 

28 Feb 1998 4.4 ± 0.6 7.1 ± 0.9 

18 Mar 1998 93.7 ± 20. 8 d 

05 Sep 1998 2.5 ± 0.3 2.2 ± 0.1 3.6 ± 0.4 3.0 ± 0.4 0.26 ± 0.21 

05 Sep 1999 a 2.1 ± 0.2 1.8 ± 0.1 2.5 ± 0.3 0.19 ± 0.24 

21 Dec 2000 2.6 ± 0.2 6.8 ± 1.0 

21 Dec 2001 10.0 ± 1.3 

21 Mar 2002 7.4 ± 1.6 

22 Dec 2003 5.2 ±1.1 
14 Jan 2004 < 5.9 
31 Mar 2004 6.2 ± 1.9 
01 Apr 2004 < 8.3 
05 Apr 2005 4.6 ± 1.7 

08 Apr 2005 8.0 ± 3.2 

13 May 2005 8.4 ± 1.7 

14 May 2005 10.4 ± 3.8 

15 May 2005 7.0 ± 1.9 
25 Mar 2006 < 9.4 

29 Mar 2006 < 9.4 
31 Mar 2006 < 7.1 
01 Apr 2006 < 8.6 

(Average) 2.0 ± 0.1 2.2 ± 0.1 2.7 ± 0.2 2.9 ± 0.4 6.9 ± 0.2 6.9 ± 0.6 C 0.21 ± 0.04 

(x 2 /dof) 3.80 4.92 3.22 11.69 0.05 1.06° 25.81 

(Probability) < 10~ 4 < 10~ 4 0.002 < 10~ 4 0.817 0.386° < 10~ 4 



Note. — Column 1: observation date; Columns 2-9: flux densities at 1.7, 5, 8.4, 15, 22, 43, 100, and 347 GHz, respectively; Column 10: spectral 
index, derived by the least-squares fit of S u oc u +a at 5-22 GHz. The last three lines: weighted average of flux densities at each frequency, reduced 
chi-squarc, and probability as statistics for the hypothesis of a constant value. 

a Alrcady published by Nagar ct al. (2001). 

b Aperture size of 3'.'8 (see Section 3.3). 

c Dcrived from only nine detected data sets. 

d Aperture size of 15'.'1. 



free. The typical system noise temperature was about 
140 K in the double-sided band. For gain calibration, 
we scanned between the target and a reference calibra- 
tor, 1150+497, every 20 or 25 minutes. The flux scales 
of the calibrator were determined from relative compar- 
isons with known flux calibrators such as Uranus. In 
general, the uncertainty in the flux scale is expected to 
be 10%-20% as a result of accumulating the differences 
of residual errors between a reference calibrator and a 
flux calibrator. These residual differences can be due to 
pointing error, seeing effect, calibrations of atmospheric 
attenuation and elevation-dependent antenna efficiency, 
and temporal variations in polarization and flux density 
of calibrators. To avoid the residual errors of calibra- 
tion from differential atmospheric attenuation and an- 
tenna efficiency, we observed 1150+497 and a flux cal- 
ibrator at almost same elevations within 1 hr. These 
observations were made within a few days of each other, 
typically, or were separated by up to 10 days in worst 
cases (from observations for NGC 4258); 1150+497 is a 
quasar with a large amplitude but gradual light curve 
on the timescale of years, which is much longer than the 



intervals between NGC 4258 and flux calibration obser- 
vations. In the present study, we assume a 15% error in 
the flux scaling factor. We observed bright quasars for 
about 30 minutes for bandpass calibration. 

The data were reduced using the UVPR0C-II pack- 
age (Tsutsumi et al. 1997) by standard means, includ- 
ing flagging bad data, baseline correction, opacity cor- 
rection, bandpass calibration, and gain calibration. Each 
daily image was synthesized from visibilities in natural 
weighting and deconvolved using IMAGR task of the As- 
tronomical Image Processing System (AIPS) developed 
at the NRAO. 



2.2. VLA Data 

In the VLA archive, we searched data sets of obser- 
vations carried out in the continuum mode using multi- 
frequency bands with high spatial resolutions of A or 
B configurations. We retrieved (1) all available data sets 
consisting of three or four bands in the range of 5-22 GHz 
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(AA1.3-6 cm) 12 for revealing spectral time evolution, (2) 
all available data sets including the 43 GHz (A7 mm) 
band for the probe of millimeter spectrum, and (3) a 
data set including a deep imaging with A-array configu- 
ration at 1.7 GHz (A18 cm) to evaluate the contamina- 
tion from extended components of "the anomalous arms" 
in the photometric measurements at higher frequencies. 
Observations of each data set were carried out at multi- 
frequency quasi-simultaneously within 3 hr typically, or 
longer (but within 1 day at most) hours for VLA data ob- 
tained as part of VLBI observations (BR0043, BG0062, 
BM0056, and BM0112). 

We reduced the data using the MPS by following the 
standard procedure. 3C 286 scans served as flux scal- 
ing factors. We followed the guidelines for accurate flux 
density bootstrapping (Per ley & Taylor 2003), includ- 
ing correcting the dependence of the gain curve on el- 
evation for the antennas and atmospheric opacity and 
by using the clean component models of 3C 286 in self- 
calibrations. The bootstrap accuracy should be l%-2% 
at 20, 6, and 3.6 cm, and perhaps 3%-5% at 2, 1.3, and 
0.7 cm (Perley & Taylor 2003). For flux scaling factors, 
we assumed uncertainties of 2% and 5% at 1.7-8.4 GHz 
and 15-43 GHz, respectively. All target images were syn- 
thesized with natural weighting and deconvolved with 
the MPS IMAGR task. 

2.3. JCMT Data 

We retrieved archival data for the 1998 March 18 ob- 
servations of NGC 4258 that were acquired by the JCMT 
at 347 GHz (A850 /im) in jiggle-map mode with the 
Submillimetre Common-User Bolometer Array (SCUBA; 
Holland et al. 1999). The data were processed with 
the SCUBA User Reduction Facility (SURF) package. 
We applied standard reduction procedures including flat 
fielding, flagging of transient spikes, correcting the ex- 
tinction, pointing correction, sky removal, and flux den- 
sity scaling. A Uranus scan gave us a flux scaling fac- 
tor, with an uncertainty assumed to be ~15% and an 
effective beam size of 15'.' 1. On the blank sky, the mea- 
sured root mean square (rms) of the image noise was 
10.5 mjy beam -1 . 

3. RESULTS 

The results of the flux measurements are presented in 
Table 2; Figure 1 shows the cm-to-submm spectrum of 
the NGC 4258 nucleus. 

3.1. Upper Limit of Dust Contamination 

The 347 GHz (A850/mi) JCMT image of the nuclear 
region with a field of view of ~ 2' shows a structure elon- 
gated along the major axis of the host galaxy with no 
apparent nuclear prominence (Figure 2). Such a struc- 
ture is very similar to that of carbon monoxide (CO) line 
emission (Regan et al. 2001). These results suggest that 
most of the submillimeter emission does not originate 
in the AGN continuum but in interstellar dust associ- 
ated with the host galaxy. We measured the flux den- 
sity from the region corresponding to the JCMT beam of 
15'.'1 (^530 pc) centered at the nucleus. The continuum 
emission may be contaminated by the CO(3-2) emission 

12 One of the data sets (AN8500) have already been published 
by Nagar et al. (2001). 
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Fig. 1. — Centimeter to submillimeter spectrum of the nuclear re- 
gion of NGC 4258. The data from VLA at 5-43 GHz (beam sizes of 
~ 0y05-3'.'8) and NMA at ~ 100 GHz (~ 3"-10") are represented 
as filled circles. Gray crosses connected with gray solid lines repre- 
sent weighted averages at each band. The dashed curve represents 
a dust continuum spectrum (see Section 3.1) determined from the 
JCMT data at 347 GHz (15'/1). The dot-dashed line represents an 
estimated upper limit of contribution from the extended emissions 
of the anomalous arms, which was derived from a flux measure- 
ment in a 3'.'8 aperture at 1.7 GHz (open square; see Section 3.3). 
The 3cr upper limits are represented by downward arrows. Upper 
limits at 115 and 225 GHz are from literature (Section 3.2) 
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Fig. 2. — JCMT image of the nuclear region (~ 4.3 X 4.3 kpc) 
of NGC 4258 at 347 GHz (A850 /an). Contour levels are crx (4, 5, 
6, 7, 8, 9), where a = 10.5 mjy beam -1 is the rms image noise. 
The cross symbol and its length represent the optical position of 
the NGC 4258 nucleus and pointing accuracy of the JCMT, respec- 
tively. The beam size (15'.' 1) is shown in the lower-left corner of 
the image. 

of only - 1 mJy in the SCUBA bandwidth of - 40 GHz, 
according to the velocity-integrated intensity peak of the 
CO(1-0) map (Heifer et al. 2003) and the intensity ratio 
CO(3-2)/CO(l-0) (Mao et al. 2010). 
In general, the dust spectrum of external galaxies can 
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Fig. 3. — VLA image with A-array configuration of the nuclear 
region of NGC 4258 at 1.7 GHz (A18 cm) in natural weighting. 
Contour levels in black curves are separated by factors of \f2 begin- 
ning at 3<r of the rms noise Ict = 42 /xjy beam -1 . The gray curves 
represent 2cr contours. The beam size of 1'.'2 X 1'.'4 at P. A. = 82?6 
is shown in the lower-left corner of the image. 

be well represented via one- or two-temperature dust 
spectral models by using the modified blackbody emis- 
sion S v oc z^-B(Td), where -B(Td) is the blackbody spec- 
trum at temperature Td, with /? ~ 1.6-2 and two tem- 
peratures of Td ~ 20 and ~ 43 K; the lower-temperature 
component dominates at < 1000 GHz in total flux (Temi 
et al. 2004, and references therein). We adopt /3 = 1 as 
an extreme case for conservative constraint at frequen- 
cies lower than 347 GHz of the JCMT data. As shown 
by a dashed line in Figure 1, the dust contribution is 
expected to be < 2.4 mJy at 100 GHz within a 15'.'1 
region, which means a strong upper limit for dust con- 
tamination toward the active nucleus of NGC 4258 in 
the flux measurements with smaller beam sizes of NMA 
(< 10") and VLA (< 3"8) at the lower frequencies (cf. 
Doi et al. 2005b, 2011). 

3.2. Nuclear Millimeter Emission 

NMA images at 100 GHz (A3mm) with a field of view 
of ~ 1' revealed only a point-like single emission compo- 
nent at the nucleus. We detected the emission for nine 
out of fifteen epochs. The flux densities were measured 
in the image domain by elliptical Gaussian profile fitting 
with the AIPS JMFIT task. We determined the total er- 
rors in the flux measurements from the root sum square 
of errors of the Gaussian fitting and flux scaling (15%; see 
Section 2.1). For detected epochs, the millimeter contin- 
uum emission is significantly stronger than the centime- 
ter emission. 

The upper limit of dust contamination is < 2.4 mJy at 
100 GHz in a ~ 15" region at the nucleus (Section 3.1). 
Moreover, an upper limit at 115 GHz was obtained on 
1998 April 20 with a smaller beam size of 6"1 x 5" 4 
(Heifer et al. 2003); an upper limit at 225 GHz was ob- 
tained on 2004 May 4 with an even smaller beam size 
of 3'/0 x 2'.'0 (Sawada-Satoh et al. 2007, see Figure 1), 
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Fig. 4. — Radio spectra for the NGC 4258 nucleus at 5-43 GHz 
obtained with the VLA A- and B-array configurations. The dashed 
line represents the fit to a power-law spectrum (SiGHz = 1-2 mJy, 
a = 0.32 ± 0.07, where S„ = S 1G iiz^ +a ) for all data from 5 to 
22 GHz. 

indicating < 0.8 mJy at 100 GHz (assuming /3 ~ 1) as 
a conceivable compact dust component. On the other 
hand, the contamination from the low-frequency anoma- 
lous arms was estimated to be < 0.5 mJy at 100 GHz 
on the basis of a spectral extrapolation using a flux den- 
sity of 3.4 mJy at 4.9 GHz for the nuclear region of 14" 
(Hummel et al. 1989) and a spectral index of —0.65 of 
the anomalous arms (Hyman et al. 2001). Consequently, 
the contamination both from dust and anomalous arms 
can be practically neglected in our NMA measurements. 
The detected millimeter emission is presumably domi- 
nated by the active nucleus. 

In a statistical sense, millimeter flux variability is un- 
clear (cf. the last line in Table 2), which may be due to 
large uncertainty in flux scaling (see Section 2.1). 

3.3. Spectral Evolution at VLA Wavelengths 

In VLA images at 5-43 GHz (AA6-0.7cm), we iden- 
tified a single emission component at the nucleus in all 
cases analyzed. No extended radio structure was appar- 
ent in these images with limited dynamic ranges. The 
flux densities were measured with the AIPS JMFIT. The 
total errors were determined from the root sum square 
of errors of the Gaussian fitting and flux scaling (Sec- 
tion 2.2). In all cases, the emissions were almost unre- 
solved; their deconvolved sizes were not beyond half of 
the FWHM of synthesized beams. The deep VLA im- 
age at 1.7 GHz (A18cm; Figure 3) exhibits compact core 
+ diffuse structures of the bases of "anomalous arms" 
to north-south directions. This low-frequency image al- 
lowed us to evaluate the contamination from extended 
components with a steep spectrum into our photometry 
toward the nucleus at higher frequencies. Photometry 
at 1.7 GHz with a 3"8-aperture (larger than the largest 
beam size in VLA observations at 5-43 GHz; see Table 1) 
resulted in 2.2 ± 0.1 mJy, which must include the emis- 
sions from both the nucleus and the bases of anomalous 
arms. The two-point spectral index between 1.7-22 GHz 
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Fig. 5.— Light curves of the NGC 4258 nucleus at 5-22 GHz in 
the first six epochs (and one more epoch at 22 GHz) with the VLA 
A-array configuration. No 22-GHz measurement was obtained at 
the fifth epoch (the 35th day: 1997 January 4). 

(1997 January 7) is inverted (a = 0.20 ± 0.11), indi- 
cating a minor contribution from the anomalous arms 
(a = —0.65; Hyman et al. 2001) in the central region of 
3'/8. Hence, 2.2 x 0/1.7 GHz) Q mJy and a = -0.65 pre- 
sumably provides very conservative upper limits for es- 
timating contaminations at higher frequencies as shown 
by a dot-dashed line in Figure 1. 

The VLA spectra for each epoch are shown in Fig- 
ure 4 and exhibit clear spectral variability with respect 
to a power-law spectrum (a — 0.32 ± 0.07) represented 
by the dashed line, which is determined by the least- 
squares ht to all epochs in the range 5-22 GHz. For 
each epoch, we also derived power-law spectral indices 
by least-squares fitting of quasi-simultaneously obtained 
data from 5 to 22 GHz (Column 10 in Table 2), re- 
sulting in from a = 0.06 to a = 0.64; the average is 
(a) = 0.21 ± 0.04. The spectral index is not stable in a 
statistical sense (a probability p < 10~ 4 for the hypothe- 
sis of a constant value; see Table 2); however, it remains 
inverted (i.e., a > 0) throughout all epochs. Two-point 
spectra including 43 GHz (a = 0.71 ± 0.54 between 22- 
43 GHz on 1998 February 28 and a = 0.59±0.30 between 
8.4-43 GHz on 2000 December 21) produce evidence of 
this trend continuing up to millimeter wavelengths. 

The spectral features and their time evolutions were 
apparently complicated. Inverted spectra and, at the 
same time, flat and marginally steep spectra were lo- 
cally observed in the frequency domain (Figure 4). The 
first six epochs show spectral features changing in only a 
month; the light curves for each frequency are shown in 
Figure 5. The flux variations of the two higher frequen- 
cies, 15 and 22 GHz, a gradual increase, appeared in the 
latter half of the six epochs. On the other hand, flux 
variations of the lower two frequencies, 5 and 8.4 GHz, 
also correlated with each other but in a different way 
from those at the higher frequencies; a modest peak-out 
is apparent around the second or third epoch, and subse- 
quently a (possible) sudden enhancement appears at the 



sixth epoch at both frequencies. We estimate the size 
of the emitting region from L < c \ dh\S v / dt\ , where 
L is the size, c is the speed of light, is the flux den- 
sity, and t is the observation time, assuming an expo- 
nential light curve (Burbidge et al. 1974; Valtaoja et al. 
1999). The synchronized increments by ~ 65% in light 
curves at 15 and 22 GHz during the last two weeks re- 
sult in L <0.02 pc (<4000 AU or < 5300Rs). A flux 
decrease at 22 GHz over 13 days from 1997 March 6 to 
March 19 (beyond the time range of Figure 5) results 
again in L < 0.02 pc. The flux density ranges and the 
timescales of variability are consistent with previous ob- 
servations at 22 GHz (Appendix A). 

4. DISCUSSION 

We found (1) a complexly variable radio spectrum and 
(2) the trend of a slightly inverted spectrum continuing 
throughout from centimeter (5 GHz, A6cm) at least up to 
millimeter (100 GHz, A3mm) regimes for the NGC 4258 
nuclear radio continuum component. 

4.1. Radio Emitting Site 

We discuss the physical origin of observed radio emis- 
sions in our data. VLA beam sizes of ~ 0'.'05-3'.'8, cor- 
responding to 1.8-130 pc, decrease with increasing fre- 
quency; the observed inverted spectra are in opposite 
sense to the resolution effect. It means the radio emis- 
sion should be dominated by a compact component sig- 
nificantly smaller than the beam sizes. VLA and VLBI 
flux densities are almost the same (e.g., ~ 3 mJy at 
22 GHz Herrnstein et al. 1998), indicating the observed 
VLA flux densities originate in a sub-pc scale region. 
On the basis of VLBI images, the radio continuum flux 
comes from the northern side of the two-sided nuclear 
jet ~ 5-10 mas away at 1.5 GHz (Cecil et al. 2000) and 
- 0.4 mas (0.015 pc or 4000Rs) away at 22 GHz (Herrn- 
stein et al. 1997) with respect to the dynamical center of 
the nearly edge-on viewed Keplerian disk of rotating wa- 
ter masers. Thus, most of the detected flux at 5-22 GHz 
in the present study is presumably of nuclear jet origin 
(rather than accretion flow), particularly from the north- 
ern side of jet, at locations > 4000Rs away from the black 
hole. 

Then, we consider the locations of 22-100 GHz radio 
emitting sites. Our finding of a slightly inverted spec- 
trum is consistent with the picture of a relativistic plasma 
flow with an opacity gradient (Blandford & Konigl 1979). 
The frequency dependence of the position of the "core," 
which is identified by an intensity peak at the root of the 
jet (i.e., the "core-shift effect"; Lobanov 1998), is fre- 
quently observed for radio-loud AGNs (Kovalev et al. 
2008; O'Sullivan & Gabuzda 2009; Hada et al. 2011; 
Sokolovsky et al. 2011) and also for the LLAGN M81 (Bi- 
etenholz et al. 2004; Marti- Vidal et al. 2011). At a given 
frequency v, the separation of the core from a central 
engine satisfies r — Vlv k with k ~ — 1 in most cases. The 
northern offset of ~ 0.4 mas at 22 GHz in the NGC 4258 
nuclear jet implies f2 s» 9 mas GHz (~ 0.3 pc GHz or 
~ 8 x 10 4 Rs GHz) with the assumption k = — 1. Ac- 
cording to this dependence, an expected position offset 
of ~ 6 mas at 1.5 GHz is actually consistent with the 
observed offset of ~ 5-10 mas at 1.5 GHz. Our finding 
of a continuous trend as an inverted spectrum at least up 
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to 100 GHz suggests that the radial profile of jet struc- 
ture maintains over the range of radio emitting sites of 
1.5-100 GHz or higher because both the frequency de- 
pendence of core shift and spectral index are relevant to 
the radial profile of jet structure (Blandford & Konigl 
1979). The emitting site at 100 GHz is, therefore, ex- 
pected to be ~ 800Rs away, which is ~ 4.5 times as 
close as the 22 GHz emitting site to the central engine; 
it is still far from an effective radius of < lOORs for an 
optically thick radio emitting photosphere of the ADAF 
(Mahadevan 1997; Herrnstein et al. 1998). 

Finally, we discuss emitting sites at even higher fre- 
quencies. We cannot determine a spectral profile at 
> 100 GHz because the JCMT image at 347 GHz indi- 
cates heavy contamination by dust emission toward the 
nucleus (Section 3.1). On the other hand, infrared flux 
measurements by isolating the nucleus with high angu- 
lar resolutions (Chary et al. 2000) are very suggestive. 
A flux density of 435 mJy at 17.9 /xm (16,800 GHz) is 
much higher than an extrapolation from the radio spec- 
trum of a « 0.3. It invokes a different origin, although 
the infrared-optical spectrum of a = —1.4 ± 0.1 indi- 
cates a nonthermal origin (Chary et al. 2000) as well. 
Yuan et al. (2002b) proposed a jet-dominated SED model 
for the observed broadband SED that peaks in the in- 
frared regime for NGC 4258. The SED model ascribes 
the infrared peak and variable X-ray spectrum to syn- 
chrotron and sclf-Comptonized emissions, respectively, 
from the jet-base with a shock diameter of 5Rs, which 
is a preacceleration region in the jet connected with the 
innermost accretion flow. Radio emissions cannot be at- 
tributed to the jet-base because of strong synchrotron 
self-absorption in a very compact region. An outer jet 
component as a postacceleration region in the jet on a 
scale of > lOORs (Yuan et al. 2002b) is required to ex- 
plain the radio emissions. Using our finding, the VLBI 
results, and the SED model, two predictions can be made 
for future observations: (1) the radio spectral variation 
will not be so promptly correlated with those in the in- 
frared or X-ray regimes (cf. Miller et al. 2010, for M81) 
because of significantly different linear scales of emitting 
sites and (2) the outer jet component will provide a con- 
tribution comparable with that of the jet base at submil- 
limeter wavelengths, which implies a spectral upswing 
from a « 0.3 at lower frequencies to a < 2.5 at higher 
frequencies rather than a spectral turnover. 

4.2. Mechanism of Radio Spectral Variation 

The complexly variable inverted spectrum we found 
should reflect the nature of the northern side of the nu- 
clear jet, where a large fraction of radio emission arises 
(Section 4.1). The time variation of bright systemic 
masers is strongly correlated with the northern jet emis- 
sion; it is naturally explained by intrinsic variability in 
the background continuum of the southern jet and only 
mildly saturated masers (Herrnstein et al. 1997). It 
means that the observed radio properties in the northern 
jet mainly originate in intrinsic ones, rather than in the 
external effect such as free-free absorption. On the ge- 
ometry of a warped disk model (see fig. 10 in Herrnstein 
et al. 2005), the northern jet is behind the disk at radii 
of > 0.29 pc in line of site, where intervening medium is 
not ascribable to an effective free-free absorber in terms 
of the opacity (Neufeld & Maloney 1995) or timescale 



(an inhomogeneity of ~ 10 15 cm in size and a rotating 
velocity of < 750 km s _1 ; Fruscione et al. 2005). 

In this framework, the frequency dependence of the 
core shift in the northern jet (Section 4.1) is ascribable to 
synchrotron self-absorption; the emitting site at 22 GHz 
would be almost opaque at 5 GHz. The light curves at 
lower two frequencies were very similar; those at higher 
two frequencies formed in synchronization. However, the 
former pair showed a different trend from the latter pair 
(see Figure 5 in Section 3.3). This behavior is consistent 
with a picture of different emitting regions at different 
frequencies. The observed complexly- variable inverted 
spectrum in the NGC 4258 nucleus can be interpreted 
as the superposition of emission spectra originating at 
different locations in the nuclear jet with frequency- 
dependent opacity. Even if a phenomenon propagates 
from upstream to downstream at the speed of light, a 
variation at 15 and 22 GHz would appear as a corre- 
sponding variation at 5 and 8.4 GHz a few months later. 
Such behavior was not detectable in our data for the 
first month (Figure 5); a series of the later epochs (1997 
March onward) were too sparse. 

4.3. Comparisons with M81 

Our findings of the radio spectral properties for the 
NGC 4258 nucleus allow us to understand that the under- 
lying mechanism in the jet component (Sections 4.1 and 
4.2) is basically identical to that of M81. This is because 
the two LLAGNs share all the following observational as- 
pects in radio regimes: (1) a variable and complex spec- 
tral profile, (2) slightly inverted spectrum with a rts 0.3 
in average, and (3) core shift on a nuclear jet (see Sec- 
tion 1 for M81). Presumably, the underlying mechanism 
must be compatible to other more distant LLAGNs (spa- 
tially unresolved so far) too, whose nuclear components 
are rapidly varying radio emissions (Anderson & Ulves- 
tad 2005) showing flat /slightly inverted spectra (Nagar 
et al. 2001; Doi et al. 2005b). 

Now, the two well-studied LLAGNs enable us to supply 
a physical understanding by means of quantitative com- 
parisons. NGC 4258 exhibits aspects that differ from 
M81 in terms of: (1) the apparent flux ratio of jet to 
counter jet Rp: two-sided morphology in NGC 4258 
(i?F ~ 6; Herrnstein et al. 1998), whereas one-sided in 
M81 (i? F > 170), (2) radio loudness (radio-to-X-ray lu- 
minosity ratio) Rx = vP la( n /Lx (Tcrashima & Wilson 
2003): a much radio-louder SED for M81 (R x ~ 2xl0~ 3 ) 
compared to NGC 4258 (Rx ~ 4 x 10~ 5 ), and (3) mag- 
nitude of core shift: ft w 9 mas GHz (<~ 0.3 pc GHz or 
- 8 x 10 4 Rs GHz) for NGC 4258 (Herrnstein et al. 1997; 
Cecil et al. 2000; see Section 4.1) is much larger than 
ft = 1.75±0.20 mas GHz (0.031 pc GHz or 4500Rs GHz) 
for M81 (Marti- Vidal et al. 2011). To calculate R F of 
M81, we adopted a peak intensity of <~ 100 mJy beam -1 
and three times the image noise of o <~ 0.2 mJy beam -1 
in the images at 8.4 GHz by Bietenholz et al. (2000). To 
calculate Rx, we used observed flux densities at 22 GHz 
of typically - 3 mJy and - 100 mJy for NGC 4258 
(Table 2) and M81 (e.g., Bietenholz et al. 2000), respec- 
tively, as representative of radio data because of similar 
spectral indices (a ~ 0.3) and 2-10 keV X-ray luminosi- 
ties typically ~lx 10 41 erg s -1 and ~ 2 x 10 40 erg s _1 
for NGC 4258 (Makishima et al. 1994; Fruscione et al. 
2005; Yamada et al. 2009) and M81 (Ishisaki et al. 1996; 
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Fig. 6. — (1) Boosting factors of Doppler beaming in flux den- 
sities for continuous jets. Solid and dashed curves represent the 
cases of inclinations of 14° and 82°, respectively (Equation (1)). 
Filled circles represent possible solutions for M81 and NGC 4258 
(Section 4.3). (2) Relationship between bulk Lorentz factors of 
M81 and NGC 4258. The solid curve represents Equation (2), a 
constraint from the ratio of radio loudness of M81 and NGC 4258. 
The dashed line represents Equation (3), a constraint from the 
ratio of observed core shifts of M81 and NGC 4258. The range 
(r(M8l) < 2.8) illustrated by dot-dashed curves would not be al- 
lowed because of the observed one-sidcdness. 

Markoffet al. 2008; Miller et al. 2010), respectively. In 
the present paper, we demonstrate that the differences 
in (l)-(3) can be attributed to inclination effects by fol- 
lowing discussions; M81 is presumably a pole-on viewed 
system at an inclination of 9 = 14° (Devereux et al. 
2003), whereas NGC 4258 at 6 = 82° is viewed nearly 
edge-on (Herrnstein et al. 1996). 

First, we examine the apparent flux ratio of jet to 
counter jet Rf- The Doppler beaming effect would make 
emission boosted by a factor of 5 2 ~ a in a continuous flow 
(Ghisellini et al. 1993), where the Doppler factor and 
bulk Lorentz factor are defined as S = \T(1 — f3 cos 9)]~ x 
and r = (1 — /3 2 )~ a5 , respectively, and /3 is the jet speed 
in the unit of speed of light. The apparent flux ratio of 
jet to counter jet is expressed as 



1 + /3 cos 9 
1 - p cos 9 



1-a 



(1) 



With this equation, 9 = 82° for NGC 4258 implies 
Rp < 1.61 at any value of T (i.e., almost symmetric); the 
observed value of Pp ~ 6 must be affected by free-free 
absorption on the southern jet (Herrnstein et al. 1996). 
On the other hand, M81 with R F > 170 requires T > 2.8 



at 9 = 14°, and then, 5 > 3.8. Thus, we can understand 
the two-sidedness in NGC 4258 as an inevitable result, 
and obtain a constraint of jet speed for M81. 

Second, we verify the difference in radio loudness Rx- 
At these inclinations, the Doppler beaming effect makes 
the (approaching) jet emission boosted for M81 and de- 
boosted for NGC 4258 at any bulk Lorentz factors (Fig- 
ure 6, panel 1). Therefore, we can naturally under- 
stand that M81 is significantly radio louder in SED than 
NGC 4258, if radio and X-ray emissions originate in a 
beamed and a non-beamed component, respectively. In- 
deed, studies of SED modeling ascribed radio and X-ray 
emissions to the outer jet (synchrotron) as a postacceler- 
ation region and the jet-base component (inverse Comp- 
ton) as a preacceleration region, respectively, in both 
cases of NGC 4258 (Yuan et al. 2002b) and M81 (Markoff 
et al. 2008). As a result, for NGC 4258, the jet-base com- 
ponent dominates throughout from IR and X-ray (Yuan 
et al. 2002b), which is presumably due to deboosting on 
the outer jet. On the contrary, the outer jet plays a major 
role in the SED of M81 (Markoff et al. 2008) . Here we at- 
tempt to constrain the bulk Lorentz factors of the nuclear 
jets (corresponding to outer jets in the SED models). We 
consider an observed radio power as a boosted emission 
(P° bs = 5 2 ~ a P™ t ) at a given frequency v. We assume 
that the intrinsic radio power (P™*) is relevant to a to- 
tal radiated synchrotron luminosity {L syn ). Moreover, 



we adopt a tentative assumption of L s 



oc 



OC Ly 



where Lj ct and Lx are the total jet luminosity and X-ray 
luminosity, respectively, as non-beamed (T ~ 1) compo- 
nent. Then, 
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Under fixed inclinations 9, Lorentz factors T 
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(2) 



and 



r(NGC4258) in Doppler factors S can constrain each other 
as represented by a solid curve in Figure 6, panel 2. We 
can find T 



(NGC4258) 



> 3 based on Equation (2) together 
with the constraint from the observed one-sidedness in 
the M81 nuclear jet (the dot-dashed curve in Figure 6, 
panel 2 is in the forbidden region). 

Third, we discuss the difference in magnitude of core 
shift il. The value of f2 also depends on both inclination 
angle 9 and Lorentz factor T. We employ Equation (11) 
in Lobanov (1998) and then consider the ratio of fl for 
NGC 4258 and M81: 

z)sm 5 / 3 9 f t x \3/3 



n oc 



i 



\tVt 



(NGC4258) 



(M81) 



10, 
9 mas GHz 
1.75 mas GHz ' 



(3) 



where Z?l is the luminosity distance. We assume that (f> 
and O, the opening angle of jet and the range of emission 
region along jet (O = In (r max /r m j n ); Blandford & Konigl 
1979), respectively, are the same between NGC 4258 and 
M81, and then eliminated in the ratio in Equation (3). 
We assume L s 



syn oc Lx again. Under fixed inclinations 9, 
Lorentz factors F( M81 ) and r( N Gc 42 5g) involved in Equa- 
tion (3) can constrain each other as represented by a 
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dashed line in Figure 6, panel 2. The dependence sug- 
gests that the jet speed is not so different but may be a 
slightly lower speed in NGC 4258. 

Thus, the combination of these three ratios in terms 
of (1) i?F, (2) Rx, and (3) fi constrains Lorentz factors: 
the intersection of the two equations (Equations (2) and 
(3)) is at r (NGC4258) - 2.8 {(3 ~ 0.93 2 and S ~ 0.42) 
and F( M81 ) <~ 3.7 ((3 ~ 0.963 and S ~ 4.1), shown as a 
filled circle in Figure 6, panel 2 and corresponding val- 
ues of 5 2 ~ a in Figure 6, panel 1. This T(M8i) value is 
consistent with T > 2.8 derived from the observed one- 
sidedness using Equation (1). Thus, jet speeds in both 
M81 and NGC 4258 are fairly relativistic, although these 
discussions are based on several assumptions. Roughly 
speaking, a lower ratio in Equation (2) leads to a larger 
r(M8i)! a smaller ratio in Equation (3) leads to a larger 
T(NGC4258); r (M81) > r (NGC 4 258 ) is always seen in a rea- 
sonable range of T m 1.5-30 and a conceivable uncer- 
tainty of ~ ±50% in the ratios. If we adopt another as- 
sumption as L syn oc i^-iokoV^BH 1 (t ne revised relation 
of the fundamental plane by Plotkin et al. 2012) instead 
of L syn oc L 2 -iokcV, r( NGC4258 ) <~ 3.2 and T( M81 ) ~ 6.8 
are obtained. 

Consequently, the observed differences between 
NGC 4258 and M81 in (l)-(3) can be naturally inter- 
preted by the Dopplcr beaming effect due to fairly rela- 
tivistic jet speeds and quite different inclinations. In pre- 
vious studies, T > 2-3 was also predicted for NGC 4258 
(Falcke & Biermann 1999; Yuan et al. 2002b) and M81 
(Falcke 1996; Markoff et al. 2008) by the model of the 
outer jet freely expanding and weakly accelerated by in- 
ternal pressure gradients, which can naturally explain the 
frequency-dependent size and spectral index of compact 
radio core (Falcke & Biermann 1995). In fact, mildly or 
fairly relativistic jet speeds have been inferred from ob- 
servations from other LLAGNs/nearby Seyfert galaxies, 
e.g., Arp 102B (> 0.45c; Fathi et al. 2011), NGC 4278 
(~ 0.76c (r w 1.5); Girolctti et al. 2005), and NGC 7674 
(~ 0.92c in projection; Middelberg et al. 2004). 

5. SUMMARY 

We investigated the spectral properties of the second 
closest LLAGN radio nucleus in NGC 4258 and discussed 
the physical origin of these properties and comparisons 
with the closest LLAGN M81. The research is summa- 
rized by the following points: 

• Radio flux variability was detected at all bands in 
the centimeter regime (5-22 GHz). In a statisti- 
cal sense, millimctric variability is unclear, which 
may be due to an insufficient number of epochs 
(43 GHz) and the large uncertainty in flux mea- 
surements (100 GHz). The VLA 1.7 GHz deep im- 
age and the JCMT 347 GHz observation provided 
the upper limits of extended synchrotron and dust 
contaminations, respectively. 

• Time-averaged flux densities showed a slightly in- 
verted spectrum (a ~ 0.3) at 5-100 GHz. A spec- 
tral turnover at a higher frequency is unclear from 
the available data. 



• Quasi-simultaneous multi frequency VLA observa- 
tions at 5, 8, 15, and 22 GHz revealed a variable 
and complicated spectral profile, but always an in- 
verted spectrum in all 10 epochs. The first six 
epochs contained intra-month variation. The light 
curves at lower two frequencies were very similar; 
those at higher two frequencies formed in synchro- 
nization. However, the former pair showed a dif- 
ferent trend from the latter pair. 

• The observed radio spectral properties can be in- 
terpreted as the superposition of emission spectra 
originating at different locations with frequency- 
dependent opacity along the nuclear jet. 

• The NGC 4258 nucleus demonstrates the follow- 
ing characteristics similar to those of M81: (1) a 
complexly variable spectral feature, (2) a slightly 
inverted spectrum with a « 0.3 in average, and 
(3) apparent core shift in the AU-scale jet. 

• We made quantitative comparisons between 
NGC 4258 and M81 in terms of jet/counter jet ra- 
tio, radio loudness, and magnitude of core shift. 
We found a solution of the combination of bulk 
Lorentz factors F( N Qrj 4258 ) <~ 2.8 and F( M81 ) <~ 3.7 
for inclinations of 82° and 14°, respectively. These 
discussions are based on several assumptions; how- 
ever, these results are consistent with relativistic 
jet speeds inferred in a few other LLAGNs in liter- 
ature. 
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APPENDIX 

PREVIOUS RADIO CONTINUUM OBSERVATIONS 

Spectrum (by quasi- simultaneous observations). Ho & Ulvestad (2001) reported a slightly resolved core of Si.4GHz = 
2.73 and Ssghz = 2.18 (mJy), a = —0.18, in matched angular resolutions of ~ 1", which were in agreement with 
previous observations at a similar resolution (Vila et al. 1990). Hyman ct al. (2001) reported the time variation 
of spectral indices ranging from —0.18 to —0.41 between 1.4-5 GHz toward the nucleus using the VLA A-array 
configuration. On the other hand, the spectral index map shows a relatively uniform distribution of the steep spectrum 
of a = — 0.65± 0.10 throughout "the anomalous arms". An inverted spectrum had been seen toward the nucleus at 
higher frequencies: S5GHZ = 2.4 and Sisghz = 3.2 (mJy), a = +0.4, in matched resolutions of ~ 3" (Turner & Ho 
1994). Nagar et al. (2001) reported a folded spectral feature (2.1, 1.9, and 2.8 mJy at 5, 8.4, and 15 GHz, respectively) 
by the VLA A-array observations at matched resolutions of 0'.'5. 

Variability. Significant flux variability has been previously suggested at 22 GHz: the continuum component in VLBI 
images seems to vary by a factor of several times (S22GHZ = 0.7-4.0 mJy) during 15 epochs in the period 1997.18- 
2001.61, although these measurements are based on a limited range of baseline lengths (Argon et al. 2007). Herrnstein 
et al. (1997) also showed flux variability of the continuum component at 22 GHz from < 2 to 6 mJy based on both 
VLA and VLBI in the period 1994.3-1996.3. The emission varied by up to 100% over the timescale of a few weeks. 
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